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ABSTRACT 

We jointly model 22 early-type gravitational lens galaxies with stellar dynam- 
ical measurements using standard CDM halo models. The sample is inhomoge- 
neous in both its mass distributions and the evolution of its stellar populations 
unless the true uncertainties are significantly larger than the reported measure- 
ment errors. In general, the individual systems cannot constrain halo models, in 
the sense that the data poorly constrains the stellar mass fraction of the halo. 
The ensemble of systems, however, strongly constrains the average stellar mass 
represented by the visible galaxies to 0.026 ± 0.006 of the halo mass if we ne- 
glect adiabatic compression, rising to 0.056 ±0.011 of the halo mass if we include 
adiabatic compression. Both estimates are significantly smaller than the global 
baryon fraction, corresponding to a star formation efficiency for early-type galax- 
ies of 10%— 30%. In the adiabatically compressed models, we find an average local 
B-band stellar mass-to-light ratio of (M/L) = (7.2 ± O.5)(M /L ) that evolves 
by d\og{M / L) / dz = — 0.72±0.08 per unit redshift. Adjusting the isotropy of the 
stellar orbits has little effect on the results. The adiabatically compressed models 
are strongly favored if we impose either local estimates of the mass-to-light ratios 
of early-type galaxies or the weak lensing measurements for the lens galaxies on 
100 kpc scales as model constraints. 

Subject headings: early-type galaxies, gravitational lensing, stellar dynamics 



1. Introduction 

The presence of dark matter in early type galaxies is clear on large scales, based on 
both weak lensing (e.g. Kleinheinrich et al. 2006, Mandelbaum et al. 2006) and X-ray (e.g. 
Humphrey et al. 2006) studies. The distribution of the dark matter and the mass fraction 
represented by the stars are less well-determined because of the difficulties in measuring early- 
type galaxy structure in the transition region between the stars and the dark matter. Stellar 
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kinematic studies of the central regions, when compared to estimates of stellar mass-to-light 
ratios, have argued either that there is little dark matter inside an effective radius (e.g. 
Gerhard et al. 2001) or that there is a substantial dark matter fraction (e.g. Padmanabhan 
et al. 2004). The significance of these differences depends on the reliability of estimating 
the stellar mass from combinations of photometry, spectroscopy and population synthesis 
models. Studies on larger scales using planetary nebulae, have found examples of galaxies 
with falling rotation curves (Romanowsky et al. 2003), while the globular clusters in one 
of these systems show a flat rotation curve (Pierce et al. 2006). Surveys of structure with 
gravitational lenses (e.g. Rusin & Kochanek 2005, Treu et al. 2006) indicate that the typical 
lens has a flat rotation curve on scales of l-2R e , but the interpretation of the scatter around 
the mean structure has been used to argue for both inhomogeneity (e.g. Treu & Koopmans 
2002a, Kochanek et al. 2006) and homogeneity (e.g. Rusin & Kochanek 2005, Koopmans et 
al. 2006) in the mass distributions. 

In this paper we reanalyze a sample of 15 lenses from the Sloan Lens ACS Survey 
(SLACS, Bolton et al. 2006) and 7 lenses from the Lens Structure and Dynamics Survey 
(LSD, Koopmans & Treu 2003) that have both mass estimates from the lens geometry and 
central velocity dispersion measurements. Koopmans et al. (2006) analyzed the sample using 
a simple, global power law model, p oc r -7 , for the mass distribution to find a mean slope of 
7 = 2.01^003 where 7 = 2 corresponds to a flat rotation curve. With a nominal scatter in 
the slope of only 0.07, Koopmans et al. (2006) argue that the halo structures appear to be 
fairly homogeneous. It is difficult, however, to relate these power law models to theoretical 
models of halos or to evaluate the significance of the scatter in the slope. Additionally, the 
models are somewhat unphysical because they allow mass distributions more compact that 
the luminous galaxy. 

Here we re analyze the SLAC S and LSD lens samples using a more physical mass model 



that combines a iHernquistl (119901 ) profile for the stars with a Navarro, Frenk & White (1996, 
NFW) model for the dark matter. By comparing the mass inside the Einstein ring of the 
galaxies (M#(< Re)) with the mass needed to produce the observed velocity dispersion, we 
can estimate the stellar mass fraction and the stellar mass-to-light ratio explicitly. We use a 
Bayesian formalism that allows us to quantitatively address the homogeneity of the sample. 
We review the data and describe our mass models and analysis techniques in §2J In §31 we 
discuss the models of the individual systems ( §3.ip . the homogeneity of the sample ( §3.2[) . 
and finally, the stellar mass fraction, the mean stellar mass-to-light ratio and the evolution 
of the mass-to- light ratios ( §3.3p . We summarize our results in §H 
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2. Data and Method 
2.1. Data 

In this paper we reanalyze the data for 15 lenses from the Sloan Lenses ACS Survey 
(SLACS) and 7 lenses from the Lens Structure and Dynamics Survey (LSD). We neglect 
two lenses, Q0957+561 and Q2237+030, from the LSD, since Q0957+561 is in a cluster and 
Q2237+030 is a barred spiral galaxy, to leave us with a sample of 22 galaxies with measured 
velocity dispersions, effective radii, Einstein radii, enclosed masses and rest frame B band 
magnitudes taken from the original analyses. F or convenience, w e summarize the data in 



Tabled], particularly since the equivalent table in lTreu et al.l (120061 ) contains ordering errors. 
For consistency we have adjusted all the data to a flat ACDM cosmological model with 
Q m = 0.3, Q A = 0.7 and H = 70 km s" 1 Myr" 1 . 



2.2. Method of Analysis 



We mod el the lenses with the two component mass model for l enses introduced by 



Keetonl (120011 ). which was also used for e arly- type galaxies in the SDSS by lPadmanabhan et al 
( 120041 ) . It consists of a lHernquistl (119901 ) model for the luminous galaxy and a Navarro, Frenk 
& White (1996, NFW) profile for the dark matter halo. By using this physically motivated 
model rather than the simple power law normally used by the LSD/SLACS studies, we both 
better connect the results to theoretical halo models, and avoid models in which the dark 
matter can be more centrally concentrated than the stars. In essence, we use the mass 
enclosed by the Einstein radius of the lens to set the virial mass M V i r , the stellar velocity 
dispersion to determine the stellar mass fraction /*, and theoretical halo models to constrain 
the halo concentration c. Finally, by comparing the derived stellar mass to the observed lu- 
minosity, we can estimate the rest frame B-band mass-to-light ratio M*/L and its evolution. 



The iHernquistl (119901 ) model used for the luminous lens galaxy is defined by 
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where the scale length r H = 0.551i? e is matched to the measured de Vaucouleurs profile 
effective radius R e and the stellar mass M* = f*M vir is related to the virial mass by the cold 
baryon/stellar mass fraction /*. The NFW profile (jNavarro. Frenk fc White! Il996l ) used to 
model the initial dark matter halo is defined by 
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where the scale length r s is related to the virial radius by c = r V i r /r s , /(c) = ln(l + c) — 
c/(l + c), and M^ m = (1 — f*)M v i r is the mass in dark matter. The average concentration 
was modeled by 
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1 + z V.8.12 x 10 12 hM Q J ' ^ 
and the individual halos have a log-normal disper sion in their concentrations of a c = 0.18 
(base 10) around the average (IBullock et al.ll200ll ). These initial models neglect the com- 
pression of the dark matter density profile by the more concentrated baryons. We esti- 
mated the change s in th e dark matter distribution using the adiabatic compres sion model of 
Blumenthal et al.l (119861 ). This approximation may exaggerate the compression (IGnedin et al. 
20041 ). so we should regard our compressed and uncompressed results as bounding the pos- 
sible effects of adiabatic compression. 

The observations provide two constraints, the mass inside the Einstein radius, and the 
stellar velocity dispersion. For any value of c and /*, we use the projected mass inside 
the Einstein radius to determine M V i r (which also determines r V i r ), then use the spherical 
Jeans equation and a constant orbital isotropy (3 to compute the velocity dispersion expected 
for the measurement aperture. The effects of seeing were modeled using a Gaussian PSF 
with the observed FWHM of the observations. Given the estimated dispersion <Ti tmo dei, the 
measured dispersion Oi and its uncertainties e a i for galaxy i, we estimate a goodness of 
fit Xi 2 (. a i) — {p%,moAd — v i) 2 leii- We model the mass-to-light ra t ios of the stars using the 



standard power 
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van Dokkum & Franx 


1996; 
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2005 


Koopmans et al. 


2006), 
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where (M*/L) is the value today and d\og{M*/ L) / dz is the rate at which it changes with 
redshift z. This in turn defines a goodness of fit x 2 ((M/L)j) with which the model fits the 
logarithm of the mass-to-light ratio (M/L)i of galaxy i, defined by the ratio of the estimated 
stellar mass (a model parameter) to the observed luminosity, given its uncertainties en = 
A(log(M/L)j = ALj/(ln(10)Lj). These two terms define a probability of fitting the velocity 
dispersion P(o$|£) = exp(— x 2 (crj)/2)/v / 27re (J j and the mass-to-light ratio P((M/L)i\£) = 
exp(— x 2 ((M/L)j)/2)/\/27reLi given the model parameters /*, c, (M*/L)o and d\og(M / L) / dz 
which we abbreviate as £. Combining the two terms, we have the probability of the model 
fitting the data Di for galaxy i 

P i {D i \t) = P{a i \OP{{M/L) i \(i). (5) 



In addition to the measurement errors listed in Table [H we should also consider sources 
of systematic errors. The essence of the method is to compare the mass inside the Einstein 
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ring M(< R e ) to a virial mass estimate from the velocity dispersion a^R/G. We can identify 
five sources of systematic errors. First, while there is little uncertainty in Me, some of the 
mass may be projected surface density from either a parent group halo to which the lens 
belongs, or from another along the line of sight. The extra density, k = S/S c in dimensionless 
units, modifies the mass inside the Einstein radius by ttk,R e 'E c , so we can think of its effects 
as a systematic error in interpreting a v of e^ = a K /2. The fu ll probability distribution 



of k is skewed to positive values (e.g. iTakada fc Hamanal 120031 ) . but we will ignore this 
problem and assume a K ±2 0.05 since the positive tail of the distribution is associated with 
detectable objects (galaxies and clusters). This systematic error also affects estimates of 
the mass-to-light ratios. Second, there are 1 — 10% uncertainties in the galaxy effective 
radius measurements which contribute uncertainties of 0.5% to 5% to our interpretation 
of the velocity dispersion. Third, the measured velocity dispersion is a Gaussian fit to 
the s pectrum, which is not id entical to the rms velocity appearing in the Jeans equation 



(e.g. iBinney fc Tremaindll987l ). The difference can be estimated from the typical Gaussian- 
Hermite coefficients \hi\ — 0.02 (IBender. Saglia fc Gerhar d 1994 ) as a fractional error in a v 
of or der v6|^4| — 0.05 in the velocity dispersion (e.g. Ivan der Marel. van Dokkum fc Franx 
20031 ). Fourth, non-sphericity, (somewhat to our surprise) leads to negligible systematic 
errors provided we use the intermediate scale length (the geometric mean of the semi-major 
and minor axes), at least in the limit of the tensor virial th eorem. It leads to large errors if any 
other scale length is used. iBarnabe &: Koopmans I (120071 ) have taken the first steps towards 
removing these two dynamical problems, although they are restricted to oblate two-integral 
models which may not be appropriate for massive elliptical galaxies. Fin ally, calibration 



errors in the velocity dispersions contribute fractional errors of order 0.03 (see iBernardi et al. 



2003al ). Combining all these contributions in quadrature, which corresponds to assuming a 
Gaussian model for each systematic error, we estimate that the typical systematic uncertainty 
to interpreting the velocity dispersions is approximately 8% with the exact value depending 
on the uncertainties in the effective radius. 

Our statistical methods are chosen so that we can understand the homogeneity of the 
lens galaxies in either their evolution or their dynamical properties and estimate their average 
properties in the presence of inhomogeneities. We will analyze the results using two Bayesian 
methods. In the first method, we will fit the data while simultaneously estimating the 
systematic errors in the velocity dispersion and the mass-to-light ratio. When combined with 
the measurement errors, these define new uncertainty estimates for the data which we will call 
the "bad" case errors in comparison to the original uncertainties (the "good" case). These 
broadened uncertainties can be representative of either true systematic uncertainties, such as 
the ones we discussed above for the dynamical measurements, or indicative of inhomogeneities 
in the structure or evolution of the galaxies. In the second method we will compare these two 
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cases using the approach outlined in lPresa (119971 ) to determine the degree to which the sample 



homogeneous or heterogeneous. In this method, we assume that there are probabilities p a and 
Pl that the galaxies have homogeneous structures or evolutionary histories in the sense that 
the scatter in the measurements is simply determined by the "good" measurement errors. 
There are then probabilities 1 — p a and 1 — pi that the galaxies are not a homogeneous 
group in either their structure or their evolution, where we characterize this by assuming 
that the uncertainties in the velocity dispersion and the mass-to-light ratio are significantly 
broadened to be the "bad" measurement errors. In essence, we are determining the relative 
probabilities of the stated measurement errors and our estimate of the true uncertainties 
from the first method. Both approaches provide uncertainties on the average properties of 
the sample that account for potential inhomogeneities, although the second method is a 
better formal approach since it can reject individual objects. 

In the first approach, we will estimate the fractional systematic errors e a and in the 
velocity dispersion and luminosity. The \ 2 expressions are modified to use uncertainties of 
e a i — > ■sj e 2 ai + el of and en — > a/ e 2 Li + e\ for the velocity dispersions and the logarithm of 
the mass-to-light ratios respectively. We assume logarithmic priors for /*, (M*/L) , and 
(d(M/L) /dz), and the theoretical prior defined by Eqn. (Ej) for the concentration q. Note 
that we are forcing all galaxies to have the same concentration, which should have no sig- 
nificant impact given the scales we are studying. The priors for the systematic errors, 
P{e c ) = 1/a/ (e^) + e 2 a o\ and P{&l) = W ( e D + e L naturally switch between uniform 
priors for systematic errors small compared to the mean square measurement errors ((e^) 
and (e 2 Li )) and logarithmic priors for large systematic errors. The resulting probability dis- 
tribution for the fractional errors is then 

P(e a ,e L \D)ocP(e a )P(e L ) f diP{i)\[P (Ak,0 P (A|ez,£) (6) 

where P (Di\e a ,£) and P (Di\eL,£) are the probability distributions modified by the addition 
of the systematic errors e a and e^. We then use these systematic error estimates to define 
the uncertainties used for the "bad" case in our second formalism. 



The second, iPressI ( 119971 ) approach properly weights all combinatoric possibilities of 
the individual systems being members of a homogeneous sample or not. Let an d 
Poi((M / L)i\£) be the probabilities of the data given the parameters for galaxy i if it is a 
member of a homogeneous group based on the measured, "good" uncertainties, and PBi(&i\0 
and PBi{{M/ L)i\^) be the probabilities if it is not and we should be using the "bad" uncer- 



taintie s based on the systematic error estimates derived from our first method. The [Press 



(119971 ) provides estimates of the relative likelihoods describing either the full sample or the 



individual systems by the "good" or "bad" data model. If we want the Bayesian probabil- 
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ity distribution for the parameters £ properly weighted over all possible group membership 
combinations, we find that 



(7) 



where 

= \prPaiiai\Z) + (1 -p CT )P B i(^|0] [p L P Gl ((M/L),|0 + (l-p L )P m {(M/L)i\0] (8) 

and where P(£) sets the prior probability distributions for the parameters. We assume a 
uniform priors for p G and p^. We obtain the probability distribution for any parameter by 
marginalizing Eqn. (j7|) over all other variables and then normalizing the total probability 
to unity. We can also estimate the probability that the sample is homogeneous in either its 
structural or evolutionary properties as 

P(p*,p L \Di) cc J dZPffilUFi (9) 

and the probability that a particular galaxy is in the dynamically homogeneous class is 

P(o"j G homogeneous | .D) = — — % — — (10) 

Ai + ijj 

where 

A l = J d£P(£) J d Pa dp LPa P Gi (a t ) \p L P Gi {{M / L)i\£) + (1 - p L )P Bt ((M/L) t \0] U^jFj aud(ll) 
Bi = J d£P(0 J d Pa d PL (l - p^Pb^IplPg^M/L)^) + (1 - p^PBidM/L^U^jFj. 

A similar set of expressions gives the probability that the galaxy is in the set of galaxies with 
a homogeneous evolutionary history. 



3. Results 

We divide our discussion of the results into three subsections. First, we present the 
results for the individual galaxies. Next, we discuss the homogeneity of the structural prop- 
erties of the galaxies. Finally, we estimate the stellar mass fraction, mass-to-light ratios and 
the rate of galaxy evolution. 
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3.1. Properties of Individual Galaxies 

Figure [I] shows contours for the goodness of fit of the models to the velocity dispersion, 
measured for each galaxy as a function of the stellar mass fraction /* and the concentration c 
once we have normalized the mass inside the Einstein radius. For these calculations, we have 
included our estimates of the systematic errors in the velocity dispersions but used the stated 
uncertainties in the luminosities. Note that the dispersion measurements cannot determine 
the halo concentrations but the goodness of fit contours always pass through the region set 
by our prior on the concentration. The permitted stellar mass fractions vary widely between 
objects. Three of the 22 objects, SDSS J0737 + 321, SDSS J1250 + 052 and PG1115 + 080, 
appear to require mass distributions that are more centrally concentrated than the stars, in 
the sense that the best fits for < 1 have y 2 > 2. This is also seen in the LSD models 



for PG1115 + 080 (jTreu fc Koopmans 2002a ). where the only models consistent with both 



the lensing constraint and the estimated velocity dispersion are more centrally concentrated 
than the stars. A fourth lens, SDSS J1627— 005, is only marginally consistent with /* ^ 1 
. Of the remaining 18 galaxies, eleven are consistent with /* = 1 (Ax 2 < 1), and seven 
are not. Four of these eleven galaxies have enormous parameter uncertainties. One problem 
for many SLACS lenses is that the scales of the velocity dispersion aperture/effective radius 
differ little from the observed Einstein radius, which limits the leverage for constraining the 
mass profile. 

Figure [2] shows the goodness of fit to the mass-to-light ratio of each galaxy given the 
best fit model for the average evolution of the sample. Most of the galaxies are consistent 
with this best fit model for the mass-to-light ratio and its evolution ( §3.3p . The mass- 
to-light ratios of the sample appear to be more uniform than the dynamical properties, 
probably for the sam e reasons that there is little scatter in the fundamental plane (see 



Bernardi et al.ll2003bl ). However, there are three 3a outliers in the sample, SDSS J1420+602, 
SDSS J1250+052 and H1543+535, all of which have very low M/L ratios compared to 
the other galaxies. Note that only one of these, SDSS J1250+052, is also an outlier in 
the dynamical fits. This is not uniq ue to our approach , since our mass-to-lig ht ratio for 



H1543+535 is comparable to that in Treu & Koopmans (2004). In Figure 6 of Treu et al. 



(120061 ). they also find significantly lower mass-to-light ratios for SDSS J1420+602 and SDSS 
J1250+052 than they do for the other SLACS members. One possible solution is that the 
lens masses are significantly mis-estimated due to contamination from a group or cluster 
halo, but only H1543+535 has a neighboring, bright galaxy and it is sufficiently distant to 
only modestly perturb the estimated mass. 
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3.2. Homogeneity 

The broad uncertainties and occasional outliers mean that it is important to have a quan- 
titative approach to determining the homogeneity of the sample and to appropriately weight 
each object when determining mean properties. This is why we introduced the Bayesian 
frameworks of §2J Fig- El shows our estimates of the fractional systematic errors from our 
first analysis method (Eqn. [6]). The best fit estimates for the fractional systematic errors in 
the velocity dispersion and luminosity are e a ~ 0.1 and cl — 0.18. The reported measure- 
ment errors lie well outside the 99.7% likelihood region. For the dynamical errors, the best fit 
systematic errors are quite consistent with our prior estimates based on simple considerations 
about the dynamical data. 

Fig. H] shows that the results for the homogeneity of the sample are very sensitive to the 
assumed uncertainties. If we simply used the stated measurement errors, then the probability 
that the sample is homogeneous in its dynamical properties (i.e. that the "good" uncertainty 
estimates are correct and the scatter is due only to measurement error) is p a < 24% and 
that it is homogeneous in mass-to-light ratio evolution is pi < 14%. Many objects have 
low probabilities of belonging to either a homogeneous dynamical subset (SDSS J1627-0053 
with p a = 0.005, SDSS J1250+0523 with p a = 0.010, SDSS J0737+321 with p a = 0.010, 
PG1115+080 with p a = 0.029, and SDSS J1420+602 with p a = 0.042) or a homogeneous 
evolutionary subset (SDSS J1250+0523 p L « 0, H1543+535 p L « 0, SDSS J0912+002 
p L = 0.001, SDSS J1420+602 p L = 0.001 and MG1549+305 p L = 0.001). Not surprisingly, 
these objects are also outliers in the individual fits from the previous section. If we include 
our estimates of the systematic errors in interpreting the dynamical measurements, then the 
probability that the sample is homogeneous in its dynamical properties rises to p a > 40%, 
but the probability of a homogeneous evolutionary population remains small at pi < 14%. 
With the inclusion of the systematic error estimates, the objects with the lowest probabilities 
of belonging to the homogeneous dynamical subset are PG1115+080 (with p a = 0.48), SDSS 
J1250+0523 (with p a = 0.49), SDSS J1627-0053 (with p a = 0.52), SDSS J0737+321 (with 
p a = 0.52), and H1417+526 (with p a = 0.55). These estimates strongly indicate that 
either the SLACS /LDS lens populations are inhomogeneous or that the measurement errors 
underestimate the true uncertainties. In a few cases, these problematic lenses show some 
evidence for disks (SDSS J1420+602, MG1549+305). 

In sum, the SLACS/LDS galaxies are homogeneous in neither their dynamical nor their 
evolutionary properties if we take the measurement errors at face value. It is likely that most 
of the problem for the dynamical measurements is that the systematic errors in interpreting 
velocity dispersions are significant and need to be included in any analysis of the dynamics 
of lenses. One of these systematic errors, surface density contributions from structures 
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other than the lens galaxy, also produces systematic errors in the mass-to-light ratio, with 
a M — g k — 0.05, but this is much too small to explain the spread in the mass-to-light ratios. 
This problem is probably caused by a combination of underestimated uncertainties in the 
luminosities and true variance in the ev olutionary history of early-type galaxies. Rusin & 
Kochanek (2005) and iTreu et al.l (120061 ) had found earlier that the lens sample was better 
fit by allowing a range for the mean redshift at which the stars formed than by assuming 
a single value, and in this analysis a range of formation redshifts would lead to non-zero 
systematic errors in the mass-to-light ratio. 



3.3. The Stellar Mass Fraction and Mass-to-Light Ratio 



We can combine the galaxies to make joint estimates of the stellar mass fraction /#, 
the mass-to-light ratio M^/L at z = and its evolution. We considered both Bayesian 
frameworks so that either the uncertainties are broadened to make the results consistent 
with all the lenses (method 1, Eqn. E]) or the outliers in the sample are properly down 
weighted in the analysis (method 2, Eqn. |7j) The two methods give similar results, so we only 
present the detailed results from the second Bayesian method. Fig. O shows the estimated 
stellar mass fraction /* for both the individual galaxies and the sample as a whole, and 
for orbital anisotropies of (3 — —1/3,0 and 1/3, where [3 = 1 — a t 2 /a r 2 is related to the 
ratio of the tangential o t and radial a r velocity dispersi ons. For isotropic, adi abatically 
compressed models, we find /* = 0.056 ± 0.011, and like iKoopmans et al.l (120061 ). we find 
that that the isotropy has little effect on the inferred mass distribution. The stellar mass 
fraction is significantly lower tha n the global baryon f raction of 0.176io'gi9 from the WMAP 



CMB anisotropy measurements (jSpergel et al.l 120061 ) . If we do not include the adiabatic 
compression of the halo, the stellar mass fraction drops to 0.026 ± 0.006, again with little 
dependence on the isotropy (3. While the uncompressed models have less dark matter in 
the central regions, the total halo mass is much larger t han in the compr e ssed models. 
To the extent that adi abatic compression occurs, but the iBlumenthal et al.l (119861 ) model 
exaggerates its degree (IGnedin et al.ll2004l ). reality is intermediate to these two extremes. 



We also fit the mass-to-light ratio as log(M*/L) = log a + bz, where a = (M*/L)o is 
the mass-to-light ratio at z = and b = d(log(M*/L))/dz is its evolution with redshift. 
Fig. El shows the likelihood contours for these two parameters for both compressed and 
uncompressed models. For the compressed isotropic models, we find a = (7.2 ± O.5)M /L 
and b = —0.7 2 ± 0. 08. This agrees with the local val ue of a = (7.3 ± 2.1)M ( t ) / from 



Gerhard et al.l (120011 ) that was used by ITreu et al.l (120061 ) . It also agrees with the ITreu et al. 
( 120061 ) estimate for the rate of the evolution b = —0.69 ± 0.08. Changing the isotropy over 
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the range f3 = —1/3, 0, 1/3 has little effect, while the model without adiabatic compression 
requires higher normalizations for the mass-to-light ratio (10.0 ± 0.3) and slightly slower 
rates of evol ution. Our analys is includes 2 lenses (PG1115+080 and H1543+535) that were 
not used by lTreu et al.l (120061 ). but excluding them from the analysis has little effect on the 
mass-to-light ratios. There are no significant changes in (M*/L)o and d\og{M*/ L) / dz if we 
neglect these two lenses. 



4. Discussion 



We reanalyzed the data from the SLACS and LSD surveys of gravitational lenses with 
velocity dispersion measurements. Our mass distribution consists of a Hernquist model 
for the luminous galaxy embedded in a theoretically constrained NFW halo model. We 
investigated the homogeneity of the sample, the stellar mass fraction the local (z = 0) 
stellar mass-to-light ratio (M ± /L) and its evolution d(log(M*/ 'L)) / 'dz. As in the earlier 
study by lKoopmans et al.l (120061 ). we found that the effects of orbital anisotropy on both the 
stellar mass fraction and the mass-to-light ratio are small. 

In most central velocity dispersion measurement provides only weak a constraint 

on halo structure in the physically interesting region. Typical limits on the mass fraction 
represented by the stars have logarithmic e rrors of order 0.5 dex. While this appears to 
contradict the conclusi ons of (for example) iKoopmans et al.l (120061 ). this is not the case. 
Koopmans et al.l (120061 ) fit mass models where p oc r -7 , and find values in the range 1.8 < 
7 < 2.3. Fig. [7] shows the expected range of this slope for our models of SDSS J0037-0942, 
where we estimated the slope by fitting the p oc r~ 7 power law model to the projected mass 
distribution over a radial baseline of R e /8 to Re that approximates the leverage in using 
stellar dynamics combined with gravitational lensing to determine halo structure. For this 
typical lens, the variations in 7 of 1.6 < 7 5, 2.06 are comparable to the system-to-system 
spread in 7 observed for the SLACS systems (IKoopmans et al.l (120061 )). Thus, the spread in 
7 observed for individual SLACS/LSD lenses is comparable to the range of values found in 
our halo models, so strong conclusions about halo structure from these system will depend 
on averages over the samples rather than the results for individual lenses. 

The critical issue for determining the sample average properties is the degree to which 
the populations are homogeneous. A heterogeneous sample cannot easily be averaged to 
determine mean properties. We find the probability of homogeneity is very sensitive to the 
uncertainties in both the velocity dispersion and the luminosity. If take the measurement 
errors at face value, there is a low probability of homogeneity in either dynamical structure 
(Po- < 20%) or evolutionary history (p L < 15%). Many lenses such as SDSS J1250+052, 
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H1543+535, SDSS J1420+601, SDSS J0912+002, MG2016+112, and MG1549+305 have 
low (< 10%) likelihoods of belonging to a homogeneous sample. The primary problem is 
probably that there are significant systematic uncertainties that must be included with the 
measurement errors. Simple considerations show that typical systematic errors in interpret- 
ing the velocity dispersions should be large compared to the measurement errors, 8% versus 
5%, and adding these estimated systematic uncertainties greatly increases the likelihood 
of dynamical homogeneity. Sources of systematic error in the mass-to-light ratio are less 
amenable to simple arguments, but should certainly include the dispersion in the average 
star formation epoch of early-type galaxies fo und in an earlier a nalyses of galaxy evolution 
with lenses by Rusin & Kochanek (2005) and iTreu et al.l (120061 ). If we simply analyze the 
data to determine the most likely systematic errors, we find that we must include fractional 
systematic errors of approximately 10% in the velocity dispersion estimates and 19% in the 
mass-to-light ratio estimates in order to make the sample consistent with the hypothesis of 
homogeneity. 

Once we account for the inhomogeneity or systematic errors in the sample, we can 
evaluate sample averages that properly account for these problems. We find that the halo 
mass fraction represented by the baryons in stars is /* = 0.056 ± 0.011 if we adiabatically 
compress the dark matter and /* = 0.026 ± 0.006 if we do not. These results are comparable 
to similar the range of es timates that relied on stellar pop ulation models to estimate the 
stellar mass. For example, iLintott. Ferreras fc Lahavl (120061 ) obtained a stellar mass f r action 
of 8% by fitting monolithic collapse models to 2000 SDSS galaxies, iHoekstra et ajj (120051 ) 
found /* = 0.065^Q;QQg using weak lensing, and Mandelbaum et al. (2006) found /* = 



0.03 



+0.02 C 
_n m / 



using weak len sing. The results in these studies depend on the assumed IMF 



0.035^004 if the initial mass fraction of 



the IHoekstra et al.l (120051 ) estimate drops to 
the stars is changed from a standard Salpeter IMF to a scaled Salpeter IMF. Our results 
probably bound the stellar mass fraction since the Blumenthal et al. (1986) model we 
used may overestimate the amount of adiabatic compression (Gnedin et al. 2004). In all 
our models, the stellar mass fraction is much smaller than thecosmological baryon mass 
fraction Vt b /Vt m = 0.176±o;o?g fr o m WMAP JSpergel et al.l l2006h . which means that the 



star formation efficiency {f^Vt m /VLb) of early-type galaxies is only 15-30%. The remaining 
baryons must remain as gas distributed on the scale of the halo or its parent (group) halo. 
This discrepa ncy appears to be a common problem for any baryon accounting for normal 
galaxies (e.g. iFukugital 12004 ) and a significant constraint on star formation efficiency 



Analysis of the evolution of early-type galaxies with redshift, whether us ing the fun 



damental plane (e.g. 
mass estimates (e.g. 
(e.g. 



J0rgensen et al. 



199 



9: 



Franx et al. 



van der Marel fc van Dokkum 



20001 ; ITreu et al.ll2006h . dynamical 



2006 



i), or gravitational lens data alone 



Rusin fc Kochanek! 120051 ) . have consistently observed a steady brightening of early- 



13 



type galaxies with look-back time, albeit with modest disagreements as to the rate. Here 
we use a hybrid method, fitting the stellar mass-to-light ratios inferred from mass models of 
gravitational lenses with stellar dynamical data, to find that (M*/L) = (7.2 ± O.5)M /L 
and G?(log(M*/ 'L)) / 'dz = —0.72 ± 0.08 for the compressed models. The m ass-to-light ra- 
tio is comparable to the local value of (M»/L)n = (7.3 ± 2.1)M Q /L from iGerhard et al. 



( 120011 ) which was adopted by iTreu et al.l (120061 ) in their analysis of this data. The mass- 
to-light rat i o evo lution rate is also close to the value dlog(M*/ Vj/jiz = 0.69 ± 0.08 from 
Treu et all (20061). and marginally larger than the estimate s by Rusin fc Kochanek! (2005 ) 
and Ivan der Marel &: van Dokkuml (120061 ) . As pointed out by iRusin &: Kochanek (b()05h . the 
differences in evolution rates are partly due to different approaches to weighting the con- 
tribution of each lens to the analysis, but at least our Bayesian approach carries out these 
weightings in an objective fashion. 

In our basic analysis we cannot distinguish between the adiabatically compressed and 
uncompressed models. In essence, we can obtain the same mass distribution either using 
a high stellar mass-to-light ratio and a more extended halo or the reverse. If we impose 
the locally estimated mass-to-light ratio as a constraint, then the adiabatically compressed 
model is favored (6 to 1). We can also use mass measurements on much larger scales to 
distinguish the two models because the total halo mass is larger in the uncompressed model. 
In particular, we can calculate the weak lensing AS and compare it to the measurement 
by Gavazzi et al. (2007, also see Mandelbaum et al. 2006) for an overlapping sample of 
SLACS lenses where they found that AS = (100 ± 30)hM Q pc~ 2 on scales of 94fr -1 kpc. 
With this constraint the adiabatically compressed models are again strongly favored (by 
1000 to 1). In general, any third constraint that is dominated by the contribution from one 
mass component will break the degeneracy and lead to constraints on the degree of adiabatic 
compression or an additional structural variable such as the inner slope of the dark matter 
density distribution. 

The sample of lenses available for such analyses will continue to grow and can include 
lenses with time delay measurements (which constrain the halo structure by measuring the 
surface density near the lensed images, Kochanek 2002) as well as those with velocity dis- 
persions. With larger samples it should be possible to explore additional correlations such 
as the scaling of the stellar mass fraction and mass-to-light ratios with halo mass and the 
dependence of the evolution rate on halo mass. In the Mandelbaum et al. (2006) and 
Padmanabhan (2004) analyses of early-type galaxies in the SDSS, the changes in the mass- 
to-light ratio with halo mass are due to an increasing dark matter fraction with mass rather 
than changes in the stellar populations, but their results depend on population synthesis 



models to corre ctly estimate the stellar masses . In a large r sample of 



tested directly. Ivan der Marel fc van Dokkuml (120061 ) and ITreu et al.l (120061 ) see some evi- 



enses, this could be 
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dence for differential evolution with mass, but significantly larger samples will be needed to 
test this given the sensitivity of even the present results to sample weighting. 

We would like to that L. Koopmans, L. Moustakas, E. Rozo and T. Treu for their 
comments and R. Mandelbaum for discussions of weak lensing. 
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Note. — z\ and z B are the lens and source redshifts, R e is the lens effective radius, Re and Me are the lens Einstein radius and Einstein mass, cr ap is the measured 
velocity dispersion. We include both the measurement errors inside the listed aperture and our estimated systematic uncertainties are in (brackets). We lacked the 
seeing FWHM for SDSS J1627-0053 and SDSS J2300+0022 and simply used the average value for the other SDSS objects. / dm (< R E ) is the projected dark matter 
fraction inside th e Einstei n radius for adiaba tical ly compressed ( "c omp r" ) or not compressed ( "n o compr") models. 
Treu et alj todVl). [Koopmans et al.l ll200rj )(2l llreu et al.l ||2003|)(3), iKoopmans fc Treul ll2003h (4) 



Tonrvl 1 199£ U5),lTreu fc Koo >inan: 
Treu fc Koopmansl 11200311 ( 



Lehar et al.Nl996h fl2l 



j ll2002all(6') IKoopmans fc Treul J200: 



Treu fc Koopmansl ll2004h (9) . iKeeton et aLNiggslK lO'). iLehar et alj J 199311 f 1 1 > 
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Fig. 1. — The goodness of fit to the velocity dispersion as a function of stellar mass fraction 
/* and concentration c. The solid curves are drawn at A% 2 = l(lcr) for the fit to the velocity 
dispersion and the dotted lines are drawn at A% 2 = 4, 9 and 16 (2, 3, and 4 a). When there 
is no region with A% 2 < 1, we label the lowest contour present. The roughly vertical pair of 
solid lines indicate the la range of concentrations given the halo mass at each point. The 
inset text identifies the object, the measured velocity dispersion and the \ 2 °f the best fit. 
These are the isotropic (f3 = 0) adiabatically compressed models that include our estimate 
of the systematic uncertainties in the stellar dynamical measurements. 
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Fig. 2. — The goodness of fit to the mean trend in the stellar mass-to-light ratio as a function 
of the stellar mass fraction and the concentration c. The solid curves are drawn at Ax 2 = 1 
(la) for the fit to the mass-to-light ratio and the dotted lines are drawn at A% 2 = 4, 9 and 
16 (2, 3, and 4 a). The inset text identifies the object. These are the isotropic (f3 = 0) 
adiabatically compressed models that include our estimates of the systematic uncertainties 
in the stellar dynamical measurements. For these figures, we use the best fit evolution model 
- the fits would improve if we included the measurement errors in the evolution model. 
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Fig. 2. — continued 
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Fig. 3. — The probability distribution of the fractional systematic errors e a and in the ve- 
locity dispersion and mass-to- light ratio. The contours encompass 68%, 95% and 99.7% of the 
probability starting from the maximimum likelihoo d solution indicate d by the triangle. The 
crosse s indicate the measurement errors from from iTreu et al.l (120061 ) and iKoopmans et al. 
(l2006h . 
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Fig. 4. — The likelihood distributions for the probability that the galaxy sample is homo- 
geneous in its dynamics (p a ) or its evolution (j>l)- The contours encompass 68%, 95% and 
99.7% of the probability. The solid contours use the measurement errors for the dynamical 
uncertainties while the dashed contours include our estimate of the systematic uncertain- 
ties in the dynamical measurements. In both cases we used an adiabatically compressed, 
isotropic (/3 = 0) model. 
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Fig. 5. — The probability distribution for the stellar mass fraction /* for tangentially 
anisotropic (ft = —0.33, top panel), isotropic (ft = 0, middle), and radially anisotropic 
(ft = 0.33, bottom) dynamical models. The thin lines in each panel show the weak con- 
straints found for the individual galaxies, and the thick solid line corresponds to the joint 
probability from combining the sample. These models are adiabatically compressed using 
the modified errors and include the fit to the mass-to-light ratios. The thick dashed line in 
the middle panel shows the effect of not including the adiabatic compressions. 
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Fig. 6. — The probability distributions for the local mass-to-light ratio (M*/L) and its 
evolution d log (M*/L)/dz in the adiabatically compressed (top) and uncompressed (bottom) 
models. The contours show the 68%, 95%, and 99.7% enclosed probability contours for the 
isotropic models. The estimated evolutio n rate is marginally inconsi stent with the estimated 
of dlog(M*/L)/dz = —0.50 ± 0.19 from iRusin fc Kochanekl (120051 ) which is shown by the 
horizontal band of solid and dashed lines. The three triangles in each panel show the effect 
of changing the isotropy on the likelihood peak, with (3 = —0.33, (3 = 0, and (3 = 0.33 as 
we move from up per left to lower right. The squares with error bars are the results from 
Treu et al.l (120061 ) for the same galaxies. 
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Fig. 7. — The range for the density slope exponent 7, where p oc r -7 , for the typical lens 
SDSS J0037-0942. We estimated 7 by fitting the projected mass distribution as a power law 
between R e /8 and Re- Note that the yariation of 7 over the p hysically interesting regime is 
comparable to the scatter observed by iKoopmans et al.l (120061 ) of 1.8 < 7 < 2.3. 



